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Abstract 

A classical problem in electromagnetics concerns the representation of the electric and magnetic fields 
in the low-frequency or static regime, where topology plays a fundamental role. For multiply connected 
conductors, at zero frequency the standard boundary conditions on the tangential components of the 
magnetic field do not uniquely determine the vector potential. We describe a (gauge-invariant) consis- 
tency condition that overcomes this non-uniqueness and resolves a longstanding difficulty in inverting 
O I the magnetic field integral equation. 

1 Introduction 
o 

We consider the problem of exterior electromagnetic scattering in the frequency domain, with particular 
attention paid to the behavior of the electric and magnetic fields in the static limit. Following standard 
practice, we write E tot (x) = E in (x) + E{x) and H tot (x) = H in (x) + H(x), where {E in ,H in } describe 
known incident electric and magnetic fields, {E,H} denote the scattered field of interest, and {E tot , H tot } 
^h denote the total fields. We write Maxwell's equations in the form 

Oh 

V x H tot = -icoeE tot (1.1) 

^H V x E tot = iujfiH tot . 

CO where e, fi are the permittivity and permeability, and we define the wavenumber by k = uo^fql. The scattered 

0^ field is assumed to satisfy the Sommerfeld-Silver-Muller radiation condition: 

lim (hx*-*e)=o ( l -) . (1.2) 

CO r^oo V re/ \r ) 

o 

^S| For a perfect conductor O, two homogeneous conditions to be enforced on T, the boundary of O, are 

^H IDE]: 

n x ,E; tot = 0, n • H tot = 0. (1.3) 
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It is also well-known that, on T, we must have 

n x H tot = J, n- E tot = p/e , (1.4) 

assuming that the scattered fields are induced by a physical surface current J and a corresponding charge p 
with 

V r -J = i^p, (1.5) 

where Vr denotes the surface divergence. 

Our primary interest here is with the classical representation of electromagnetic fields in terms of the 
vector and scalar potential (in the Lorenz gauge): 

E = iuA-V(j) (1.6) 

H = -V x A (1.7) 

/i 

where A[J](x) = p f r g k (x - y)J(y)dA y and </>(x) = ^V-A, with g k (x) = £r^j. In scattering problems, 
J, the unknown surface current, is a tangential vector field. 



Using the first condition in (1.4), with care in taking the limit as x approaches the boundary T from the 



where 



exterior domain, we obtain the magnetic field integral equation (MFIE): 

i J{x) - M k [J](x) = n(x) x H in (x) (x e T), (1.8) 

M k [J](x) = n(x) x V x / g k (x - y)J{y)dA y . (1.9) 

It is a second-kind Fredholm equation, originally suggested by Maue [3] . See [4j [5j [6) [7] for detailed discus- 
sions. 

While the MFIE has a sequence of spurious resonances at higher and higher frequencies, below the 
first such resonance it yields an invertible and well-conditioned linear system, so long as the scatterer is 
simply-connected. When the scatterer is topologically non-trivial, however, the MFIE breaks down, having 
a non-trivial null-space at zero frequency. In the static limit, the dimension of the null-space is equal to the 
genus of the surface. While this problem has been carefully analyzed [HI HJ [9], no simple remedy has been 
provided to date. 

In this paper, we show that the difficulties encountered by the MFIE can be resolved through the en- 
forcement of an apparently new consistency condition on the vector potential, involving line integrals around 
B- cycles (see Fig. fll). Like the Aharonov-Bohm effect, it illustrates the fundamental role of the vector 
potential, but in a classical regime. We begin with the static case, since it is mathematically simpler and of 
importance in its own right. 

2 Magnetostatics 

In the zero frequency limit, Maxwell's equations are generally said to uncouple, with scattered electrostatic 
and magnetostatic fields denoted by (E , H ), respectively. In the unbounded domain R 3 \ Q exterior to Q, 
2£o satisfies 

Vxft = 0, V-#o = 0, (2.1) 

subject to the boundary condition 

nxE = -nxE^\ r , (2.2) 

where n denotes the outward normal to O. In R 3 \ O, H satisfies 

V x H = 0, V • H = 0, (2.3) 

subject to the boundary condition 

n.H = -n.H™\ r . (2.4) 



Spanning surfaces £2 and Si 



Figure 1: A surface of genus 2. The A-cycles are loops whose spanning surface is in the interior of the 
scatterer. The ^-cycles have spanning surfaces Si, 62 which lie in the exterior domain. The loops Li, L2 lie 
in the interior domain and pass through the A-cycles. 

Vector fields that are both curl- free and divergence- free are called harmonic vector fields. 

A fundamental difficulty arises in the static theory that is topological in nature. More precisely, if the 
number of connected components of ft is denoted by ra, then there is an m-dimensional space of nontrivial 
vector fields Eq in R 3 \ O, satisfying 

VxEj = 0, V • E% = 0, nxE% = 0\ r . (2.5) 

They are generally referred to as (exterior) Dirichlet fields and correspond to setting the electrostatic potential 
on each disjoint conductor to a different constant. This space of fields is, of course, essential in studying 
capacitance problems in a system of disjoint conductors, where Eq = — V0. 
Here, we concentrate on magnetostatics, where nontrivial solutions H (x) to 

V x H = 0, V • H = 0, n • H = 0| r , (2.6) 

are referred to as interior or exterior Neumann fields, depending on whether x G Oora; Gl 3 \0, respectively. 
We let i7 + denote the space of exterior Neumann fields and H~ denote the space of interior Neumann fields. 
It is a classical fact that the dimension of the spaces i7 + and H~ is g, where g is the genus of the boundary 

rg]. 



In short, the boundary condition (2.4) alone determines a unique field only in the simply connected 
case (g = 0). A natural representation in that setting is to seek Ho as the gradient of the magnetic scalar 
potential, H = V\I/, with \I/ a single- valued harmonic function satisfying 

n • W = -n-H [ n . (2.7) 

In the multiply connected case, additional data is needed to make the magnetostatic problem well-posed, 
all of which are designed to account for current loops that may be flowing through the handles of the domain. 
One such condition is to require that the line integrals of H around each A-cycle be specified: 



I 



H -ds = ai. (2.8) 

a. 



It can be shown that the solution to fl2.3|2.4|2.8| ) is unique (see, for example, [TOl HT1 l4l fl2l fT3l fl4] ). 



The use of the scalar potential can be extended to the multiply connected case in two ways - either 
by introducing g cuts on the boundary T and allowing for a potential jump across each cut [16], or by 
introducing g current loops in the interior domain £} that span the fundamental group of Q (essentially one 
passing through each A-cycle - see Fig. HI JTUJ [12]- In the latter case, it is straightforward to see that one 



can represent H as 

1 9 

H = V$ + -Vx^a i A[£ i ], 

^ *=i 

where Ci is a loop of unit current density flowing through Li (Fig. [I]), allowing the scalar potential \I> to be 
single- valued. 

It is also possible to extend the scalar potential approach to the full Maxwell equations at non-zero 
frequencies through a generalization of the Lorenz-Debye-Mie formalism given in [13 [18], but this involves 



non-physical variables. Unlike the above formulation, the MFIE (1.8) uses a physical unknown, the surface 



current J, and extends naturally away from the static limit through the representation (1.6)-(1.7). In this 



paper we seek to better understand, in the static case, the range and null-space of the MFIE. 

3 Magnetostatics using the vector potential 

There is a substantial literature on magnetostatics and the representation of harmonic vector fields in the 
form of a curl (see, for example, [51 [9l H2J HU IHJ EH)- We do not seek to review the theory here, except 
where it is of direct relevance to the MFIE. At zero frequency, the MFIE takes the form 



\l-M 



J = nxH in , (3.i; 



where I denotes the identity operator and we have dropped the argument x for the sake of clarity. We let 
r^~ : j = 1, . . . ,#} denote a basis for the exterior Neumann fields and {HJ 



{Hf : j = 1, . . . , g} denote a basis for the exterior Neumann fields and {H- : j = 1, . . . , g} denote a basis 



for the interior Neumann fields. We let {Z^ — Ht\ r '• j = 1, • • • ,#} denote the boundary values of the 
exterior Neumann fields and {Z~ = H~\ : j = 1, . . . ,#} denote the boundary values of interior Neumann 
fields, which are perforce vector fields tangent to T. As shown in [HI El EL the relevant null-spaces are known 
to be 

N(h-M )={nxZ+: j = l,...,g} 

)l \ (3 ' 2) 

N(^-I-M>)={Z-: j = l,...,g}, 

where M' is the adjoint of M. From Fredholm theory, the solvability conditions for the MFIE are, therefore, 
that the inner products (n x H in , Z~) = for j = 1, . . . ,#. This is always the case for an incoming field 
that is generated by currents lying outside of a simply connected neighborhood of ft [9 . In short, for the 
physical scattering problem, the operator \l — M is rank deficient but not range deficient, and it remains 
only to develop a set of constraints that make it invertible. 

We turn now to the main point of the present paper: that such constraints can be found and that they 
come from an analysis of the electric field in the limit uo — » 0. We begin by noting that on a perfect conductor, 



the vanishing of the total tangential electric field (1.3) and Stokes' theorem allows us to write 



/ E-ds = - j E in -ds = - j iufiH in -ndA 



where Sj is a spanning surface for the 5-cycle Bj. Dividing both sides by iu, using the representation (1.6), 
and noting that the line integral of a gradient vanishes, we have 

Theorem 1. Let O be a multiply- connected perfect conductor and let Bj be a B -cycle and Sj be a spanning 
surface for Bj . Then 

f A-ds = -n [ H in ndA . (3.3) 

This condition makes sense for any u > 0. Alternatively, if we have access to the vector potential A m 
which defines the incoming field, the same analysis yields the even simpler condition: 

/ A-ds = - j A in -ds. (3.4) 

Jb, Jb, 



Without Flux Condition 



With Flux Condition (21) 
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Figure 2: The induced surface current on a torus at wavenumber k = 10 -16 without (left) and with (right) 
the auxiliary consistency (flux) condition. The real part of the current in the azimuthal direction is shown. 
The incident field is due to a unit strength current loop of radius .5 located at (3, 3, 4). 
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With Flux Condition (21) 
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Figure 3: The real part of H tot ■ z on surfaces spanning the holes of the tori in Fig. [2] Note the different 
scales on the color bars, as well as the change in sign of the point- wise flux on the right. The net flux through 
the left surface is 0.019 while, on the right, it is zero to machine precision. 



Theorem 2. Let O be a multiply- connected perfect conductor of genus g with B -cycles Bi, . . . ,B g . Then 



the MFIE augmented by conditions of the form f3.3y or (3.4) for j = 1, . . . , g has a unique solution. 



Proof. As noted above, the equation (|/ — M )J — n x Hfi 1 is solvable for any physically meaningful 
right hand side. Thus the only issue is that of uniqueness; we need to show that the combined null- 



space of the operator (\l — M ) and (3.3) is trivial. The null-space of the integral operator is spanned by 



{n x Zf : j = 1, . . . ,g}. Hence, if J lies in the null-space of this integral operator, then the exterior field 



j 
takes the form 

9 



H 



^VxA[J] = £/%#+, (3.5) 



" 



see [9] . Under the correspondence between vector fields and 2-forms 

h±d Xl + /i2<9x 2 + ^3^x3 ^ h\dx2 A dx^ + /i2<^3 A dx\ + h^dx\ A dx2, (3-6) 

the fields {Hj~ : j = 1, . . . , g} are a basis for the relative cohomology group i7^ R (lR 3 \ O, T). The functionals 

H ^ j H -ndA= - j A[J]-ds for j = 1, . . . , g (3.7) 

span the dual space to this vector space. Hence if these integrals all vanish, then the coefficients {f3j} must 
all vanish as well. □ 



Remark 1. Let us suppose now that we have discretized the MFIE (1.8) with 2N unknowns (2 degrees 
of freedom for the surface current J at each of N points), resulting in the linear system Aj — b and the 
constraints \3.$ by Cj = f , where C is a g x 2N matrix. It is straightforward to show that the system 

[A + QC}j = b + Qf 

has the same solution as the constrained equation. It is invertible so long as the range of Q, a 2N x g matrix, 
has a full rank projection onto the null vectors of the adjoint A! . In the low frequency regime, it is sufficient 
to use Q — C . 

4 Time harmonic electromagnetics using the vector potential 



As soon as cj ^ 0, the MFIE ( |1.8[ ) is formally invertible and there is no need for the incorporation of 
consistency constraints. Because there is a nearby singular problem, however, the linear system is extremely 
ill-conditioned at low frequency. Without the additional constraints the condition number is 0(u~ 2 ). In that 
regime, the incorporation of the constraints improves the condition number considerably. 

5 Numerical Results 

For illustration, we consider the problem of scattering from a torus: a genus one surface of revolution (Fig. 
[2]), driven by a known current loop in the exterior. We have implemented a solver for the MFIE with and 
without the consistency condition (with a detailed discussion of the method to be reported at a later date). 
One convenient measure of the error in the solution is the flux of H tot through the hole in the torus, which 
should vanish on a perfect conductor (see Fig. [3]). 

6 Conclusions 

In this paper, we have derived a simple consistency condition (Theorem 1) for the vector potential in the 
context of scattering from perfect conductors, valid at any frequency. It is of physical interest for three 
reasons: 1) it describes an interesting correlation between the electric and magnetic field that persists at 
zero frequency, 2) it enforces uniqueness for the magnetic field integral equation (MFIE) in multiply connected 



domains in the static limit, and 3) it improves the stability and robustness of the MFIE in the low frequency 
regime. 

It is, perhaps, worth noting that the role of the vector potential in the consistency condition is, in some 
sense, dual to its role in the Aharonov-Bohm effect. In the latter case, it is the line integral of the vector 
potential around A-cycles that is critical. That integral measures the solenoidal (poloidat) current flow on 
the torus, which induces a zero electromagnetic field in the exterior. The line integral of A around ^-cycles, 
on the other hand, is sensitive to the toroidal current flow on the surface. On a perfect conductor, that 
induced current exactly cancels the flux of the incoming magnetic field through the "holes." 

7 Acknowledgments 

We wish to thank Eric Michielssen for several useful discussions. 

References 

[I] J. D. Jackson Classical Electrodynamics, (Wiley, New York, 1975). 

[2] C. H. Papas, Theory of Electromagnetic Wave Propagation, (Dover, New York, 1988). 

[3] A.W. Maue, Z. Phys. 126, 610 (1949). 

[4] D. Colton and R. Kress, Inverse Acoustic and Electromagnetic Scattering Theory, Springer Verlag, Berlin, 
1992). 

[5] W. C. Chew, Waves and Fields in Inhomogeneous Media, (IEEE Press, New York, 1995). 

[6] J. Jin, Theory and Computation of Electromagnetic Fields, (IEEE Press, Hoboken, 2010). 

[7] J.-C. Nedelec, Acoustic and Electromagnetic Equations, (Springer, New York, 2001). 

[8] K. Cools, F. Andriulli, F. Olyslager, and E. Michielssen, IEEE Trans. Antennas Propag. 57, 3205 (2009). 

[9] P. Werner, J. Math. Anal. Appl., 14, 445 (1966). 

[10] J. Van Bladel Electromagnetic Fields, (McGraw-Hill, New York, 1964). 

[II] R. Kress, Math. Meth. Appl. Sci. 8, 77 (1986). 
[12] A. Bossavit, IEE Proc. A 135, 179 (1988). 

[13] D. Mayergoyz, M. V. K. Chari, and J. D'Angelo, IEEE Trans. Magn. 23, 3889 (1987). 

[14] F. Henrotte and K. Hameyer, IEEE Trans. Magn. 39, 1167 (2003). 

[15] A. A. Blank, K. O. Friedrichs, and H. Grad, Atomic Energy Commission Report NYO-6486 (Courant 
Institute, NYU, 1957). 

[16] A. Vourdas and K. J. Binns, IEE Proc. A 136, 49 (1989). 

[17] C. L. Epstein and L. Greengard, Comm. Pure Appl. Math. 63, 413, (2009). 

[18] C. L. Epstein, L. Greengard, and M. O'Neil, Comm. Pure Appl. Math., to appear. 



